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Abstract 

IceCube, a cubic kilometer neutrino telescope will be capable of probing neutrino-nucleon in- 
teractions in the ultrahigh energy regime, far beyond the energies reached by colliders. In this 
article we introduce a new observable that combines several advantages: it only makes use of the 

G. 

upward-going neutrino flux, so that the Earth filters the atmospheric muons, and it is only weakly 
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' dependent on the initial astrophysical flux uncertainties. 
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I. INTRODUCTION 



High precision or high energy are necessary to study the matter at short distances. The 
most powerful accelerators are the cosmic ones in the outer space. We know that the Earth 
is hit by cosmic rays of very high energy, which means that there have to be astrophysical 
mechanisms capable of accelerating protons to those high energies. It is then also possible 
that the same mechanisms could produce neutrinos of high energy or that the protons 
and radiation produced interact with matter to originate extremely high energy neutrinos. 
Some candidate neutrino sources are Active Galactic Nuclei (AGN) which are the central 
regions of certain galaxies where the radiation emitted is comparable to the total radiation 
from the entire galaxy, and Gamma Ray Bursts (GRB) [2j, that are the most powerful 
explosions since the Big Bang resulting usually from the core collapse of massive stars. 

The integrated effect over all astrophysical sources in the sky where such producing 
mechanisms may operate, is expected to lead to a diffuse neutrino flux that could be detected 
by IceCube jsjj. This next-generation neutrino telescope is planned to have a good directional 
resolution, a fact that will be useful for our purposes. Another source of neutrinos that 
contributes to a diffuse flux is given by the collisions of cosmic rays with the nucleons of 
the atmosphere, however for energies above 10 5 GeV, the extraterrestrial diffuse flux should 
start to dominate over the atmospheric spectrum. 

These high energy neutrinos coming from different sources can be used to look for new 
physics effects in neutrino-nucleon interactions using the nucleons of the Earth as targets. 
In order to bound such effects, the different observables that have been studied, basically 
arise from comparing the upward-going flux that survives after passing through the Earth 
I dependent on the neutrino-nucleon cross section) to the standard model 
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prediction {4, J 



II. OBSERVING NEW PHYSICS 



In this work we define the observable a(E) in the following way. Considering only upward- 
going neutrinos, that is, the ones with arrival directions 9 such that < 9 < tt/2, we denote 
by a the angle such that the number of events for < 9 < a equals the number of events for 
a < 9 < 7r/2. Clearly, the value of a is energy dependent. At low energies the cross section 
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is lower and the Earth is essentially transparent to neutrinos. In this case, a — > 7r/3 since for 
a diffuse isotropic flux this angle divides the hemisphere into two sectors with the same solid 
angle. Obviously, for extremely high energies where most neutrinos are absorbed, a — > tt/2, 
and for intermediate energies a varies accordingly between these two limiting behaviors. 

We will consider only the diffuse neutrino flux from extraterrestrial origin and assume 
that it is isotropic. The use of the observable a(E) reduces the effects of the experimental 
systematics and initial flux dependence. The functional form of a(E) sharply depends on the 
interaction cross section neutrino-nucleon. In this conditions, if physics beyond the standard 
model operates at these high energies, it will become manifest directly on the function a(E). 

If we take into account a diffuse neutrino flux and for simplicity we consider it as de- 

n 

creasing with energy (an example is the given by the Waxman and Bacall [8J), then in a 
first analysis we can safely neglect the regeneration effects and approximate the surviving 
neutrino flux by [9|], 



where t{9) is the number of nucleons per unit area in the neutrino path through the Earth, 

i>2Re cos 8 

t{6) = N A / p{z)dz. (2) 

Jo 

Here $>o(E) is the initial neutrino flux, N A is the Avogradro number, Re is the radius of 
the Earth, and 9 is the nadir angle taken from the downward-going normal to the neutrino 
telescope. 

Upward-going neutrinos can originate, though CC v^N interactions, energetic muons that 
will traverse the detector. This is the traditional observation mode, in which the background 
due to atmospheric muons is eliminated. Simulations based on AMANDA data indicate that 
the direction of muons can be determined to sub-degree accuracy and their energy can be 
measured to better than 30% in the logarithm of the energy. The important advantage of 
this mode is the angular sub-degree resolution. On the other hand, as it was recognized 
in Ref. [a], in IceCube we will have sufficient energy resolution to separately assign the 
energy fractions in the muon track and the hadronic shower allowing the determination of 
the inelasticity distribution and the neutrino energy. Recently, the possibility to measure 
the inelasticity distribution was used to study the possibility to place bounds to new effects 
coming from leptoquarks or Black- Hole production over kinematic regions never tested before 
0, Q. In our particular case, the possibility of independently measuring the muon energy 
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and the hadronic shower energy will allow us to have a reasonable z/^-energy determination. 
Hence, in the following we shall take the z/^-energy bin partition interval as A log 10 E = 0.5. 

To a good approximation, the expected number of events at IceCube in the energy interval 
AE and in the angular interval A8 is given by 

A/" = n T T I [ d6dE u a(E v )$(E,9), (3) 

JA8 JAE 

where is the number of target nucleons in the effective volume, T is the running time, 
and (j(E v ) is the neutrino-nucleon cross section. In our analysis we are interested only in 
CC contained events, for which an accurate measurement of the inelasticity can be obtained. 
We take as the detection volume for contained events the instrumented volume for IceCube 
which is roughly 1 km 3 and corresponds to % ~ 6 x 10 38 . 

Since the definition of a is the equality between two number of events, then to a good 
approximation for each energy bin all the previous factors cancel except the integrated fluxes 
at each side. Thus, a can be defined by the equation 



/ d6sm6e-° SM( - E ^ = d9 sin Oe'^^^, (4) 

J0 Ja SM (E) 



which is numerically solved to give the results shown in 



'ig. [TJ There we have considered 
the standard model cross section as it was calculated in [ll|, and for r{9) we use Eq. [2] with 
the Earth density as given by the PREM [121 ] . 

Both the SM predictions for the total cross section and the Earth density have uncer- 
tainties that propagate into the observable a(E). In Fig. [3j where the new physics effects 
on a(E) are shown, the effects of the mentioned uncertainties of the observable are also 
included. 

To reduce background, one looks for upward-moving muons produced when neutrinos 
coming from the opposite side of the Earth interact with the nucleons in their path. Di- 
rectional reconstruction of these tracks suppresses the atmospheric muons background. By 
selecting the neutrino events that are upward-going it is possible to eliminate the atmo- 
spheric muons since they are downward-going. 
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FIG. 1: Standard Model prediction for osm(E) 



III. FOUR-FERMION INTERACTIONS 



In order to model in a general way the new physics effects on a(E), we consider general 
4-fermion interactions as given by an effective operator that includes also the SM fields 
involved in the neutrino-nucleon scattering with left-handed neutrinos. If there are new 
interactions between quarks and leptons, then the new effects should appear at an energy 
enough high. We call this characteristic energy scale for the new interactions A. At energies 
below A, these interactions are suppressed by an inverse power of A. Thus, the dominant 
effects should come from the lowest dimensional interactions with 4-fermions 131 ]. 



c = c sm + i^^^^ q^PLQi + v liJb P L v q lltl P Lqi 



(5) 



for left-handed neutrinos, where we take g 2 N = Air, and the coefficients r\ LL and r] LR can take 
the values —1, and 1. Choosing different values of A, t\ll, and tjlr, we can test the a 
observable on different scenarios of new physics. 

Using the effective operator we can calculate their contribution to the neutrino-nucleon 
inclusive cross section 



uN — ► fM + anything, 



(6) 
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FIG. 2: Feynman diagrams contributing to charged current and neutral current processes. 



where N = — - — for an isoscalar nucleon. The corresponding processes are pictured 
in Fig. [2] for charged and neutral currents. The calculation is standard and we use it to 
compare with the results of [111 ] . For the charged current, the scattering amplitude is 



where 



M 



cc 



2(Qi + M w) 



-hyPiM qrf{g' L P L + g R P R )qi 



(7) 



' , (Q 2 + M 2 W ) g 2 N 
9l = ± VLLy 



9r 



(Q 2 + My/) gl 



A 2 ~ VLR g 2 



(8) 



include the new physics effects. 
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The differential cross section for charged currents reads 



^^(w^fe)) 2 ^ 2wCC + (1 -^ CC) (9) 

+ g m cc + (i - v) 2 q cc ) 



where for an isoscalar target we have the quark distribution functions 



n cc/ n 2N _ u v (x,Q 2 ) + d v (x,Q 2 ) u s (x,Q 2 ) + d s (x,Q 2 ) 

LJ [X, LJ ) — ^ + 2 



+ s s (x,Q 2 ) + b s (x,Q 2 ) 



9l - 9l ^2 c wVll g 2 

>i i Q 2 + M 2 2 g 2 LR 
9r — 9r ^2 c wVlr-^- 



(10) 



Q cc (x, Q 2 ) = u ^Q 2 ) + ds(x ' Q2) + c s (x, Q 2 ) + t s (x, Q 2 ) 
Similarly, the neutral current amplitude is 

^ NC = -rfTMl) P7 " Pl " ^ {3lPl + g ' RPR)q " (U) 

where cw = cos 9w and 



(12) 



include the new physics effects. Here, gl = l/2 — 2/3x w ,g® = — l/2+l/3x w , g^ = —2/3x w , 
and g R = l/3xw 

The neutral current differential cross section is then 

S-^(^) , | 5,*w^-.W) (i3) 

+g£(Q l + (i-y) 2 Q% 

where the corresponding parton distributions for a isoscalar target read 
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Uf 2 _ u v (x,Q 2 ) + d v (x,Q 2 ) u s (x,Q 2 ) + d s (x,Q 2 ) 

+c s (x,Q 2 ) +t s (x,Q 2 ) 

D ^ _ u v (x,Q 2 ) + d v (x,Q 2 ) u s (x,Q 2 ) + d s (x,Q 2 ) 
v \ x i v ) — 2 2 

+ Ss (x,g 2 ) + 6 s (x,g 2 ) 

Q^, Q 2 ) = u ^Q 2 ) + ds(X,Q2) + c s (x, Q 2 ) + t s (x, Q 2 ) 
QD{ ^ g2) = u s (x,Q 2 ) + d s (x,Q 2 ) + g2) + g2) ^ 



IV. RESULTS 

In order to evaluate the impact of the observable a to bound new physics effects, we have 
estimated the corresponding uncertainties. For the number of events we have considered it 
as distributed according to a Poisson distribution and we have propagated onto the angle 
a. The number of events N as a function of «sm is 

N = 2rm T TAEa(E)(P (E) / c^sin&r^^ (14) 

Jo 

where we have considered the effective volume for contained events for which a accurate 
and simultaneous determination of the muon energy and shower energy is possible. For 
IceCube it corresponds to the instrumented volume, roughly 1 km 3 , that implies a number of 
target nucleons — 6 x 10 38 . For T we have taken an integration time of 15 yr corresponding 
to the lifetime of the experiment. 

To propagate the error on N to obtain the one on a, we note that 

AN = ^Aa, (15) 
da 

and dividing by N we obtain for Aa, 

[AN 



q sm(£) . 

a T (E)[ T (a SM (E))-r(e)] 



r #>( —, ,) 

Jo \smasM{E)J 



Aa = 

where for Poisson distributed events we have 



e 



V N 



(16) 



AN = VN. 
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(17) 



In the evaluation of the errors on a it is necessary to consider the initial flux O M to 
estimate the events rates. 

The usual benchmark here is the so-called Waxman-Bahcall (WB) flux for each flavor, 



E 



2 x l(T 8 GeV cm-V^sr -3 



(18) 



which is derived assuming that neutrinos come from transparent cosmic ray sources 
and that there is adequate transfer of energy to pions following pp collisions. However, one 
should keep in mind that if there are in fact hidden sources which are opaque to ultra-high 
energy cosmic rays, then the expected neutrino flux will be higher. 

On the other hand, we have the experimental bound set by AMANDA. A summary of 



the bounds can be found in Refs. 
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151 ] and a representative value for it is 



e I$ V am - 2 x l(T 7 GeV cm-W- 1 . (19) 

However, with the intention to estimate the number of events we have considered an 
intermediate flux level slightly below the present AMANDA experimental bound, 

E 2 U ^% ^ l(T 7 GeV cm-W- 1 , (20) 

which is the flux that we have used to estimate the uncertainties on the angle a. 

In Fig. [3] we show our results for the a observable for the most representative sets of 
parameters and the standard model prediction, including the theoretical uncertainties from 
the SM cross section, the Earth density, and the errors coming from the statistical uncer- 
tainties in the number of events. For the new physics effects we have considered the sets of 
parameters shown in Table [B 

In Fig. H] we show the differences between the values of a for different sets of parameters 
and the standard model values as a function of the energy. It can be seen that the maximum 
sensitivity is reached in the intermediate energy range (10 5 GeV < E < 10 7 GeV). In this 
figure we also include the uncertainties in the Standard Model prediction coming from the 
ones in the event numbers (Shaded region). 

In the same context, we can define another observable related to «sm(-E')- We consider 
the hemisphere < 6 < n/2 divided into two regions by the angle osm, for < 6> < ag M 
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FIG. 3: a(E) for different sets of parameters: (a) setl, (b) set 2, (e) set 5, and the standard 
model prediction with the theoretical uncertainties included (Shaded region) and the errors in a 
as resulting from the ones in the events numbers. 

TABLE I: Sets of parameters for the new four- fermion contact interactions. 



Set till Vlr A (TeV) 

111 1 

2- 1-1 1 

3- 1-1 2 

4 11 0.8 

5 -1 -1 0.8 



and 1Z 2 for ct SM < 6 < ir/2. We then calculate the ratio x between the number of events for 
each region, 

* = < 21 > 

where Ni is the number of events in the region IZi and N 2 is the number of events in the 
region 7Z 2 - By using x, the effects of experimental systematics and initial flux dependence 
are also reduced. 

If there is only standard model physics, then we have that the ratio x — 1- The new 
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FIG. 4: Differences between a for different sets of parameters and the standard model prediction 
(A a = «sm — «Setj)- We have included the uncertainties on A a calculated for the intermediate 
initial flux 0jjv (Shaded region). The bin partition energy interval is taken as Alog 10 E = 0.5 

physics effects produce a deviation from the standard value as we show in Fig. [5] for the 
different sets of parameters (different values of A, t)ll, and tjlr)- 

V. CONCLUSIONS 

In the present work we have studied a new observable that combines the surviving neu- 
trino flux after passing through the Earth in a way that reduces the experimental systematics 
and the dependence with the initial flux. This observable, the angle a, divides the upward- 
going hemisphere (with respect the arrival neutrino directions) into two homo-event sectors 
and it is dependent, of course, on the neutrino energy. The function a(E) is sharply de- 
pendent on the neutrino-nucleon cross section, which makes it a useful observable to bound 
new physics. In order to test the sensitivity of a, we have calculated the new physics effects 
coming from four-fermion contact interactions. We have also studied as another observable, 
the ratio between the number of events in the regions defined by asu{E). 

We note that the introduced observable present no deviation from the standard model 
prediction for low energies (< 10 3 GeV) at which almost no interactions occur with or 
without new physics. At high energies (> 10 9 GeV) the neutrino mean free path is so small 
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FIG. 5: Observable x{E) for different sets of parameters. 

that the integrated surviving flux of Earth skimming neutrinos equals the one corresponding 
to almost the whole hemisphere. The corresponding increase in the cross section implies that, 
given the exponential behavior of the integrated flux, a great attenuation takes place both 
taking into account new physics effects or not. 

Finally, we point out that this technique can be applied to any specific case of physics 
beyond the standard model, which is left for future work. 



Acknowledgments 

We thank CONICET (Argentina) and Universidad Nacional de Mar del Plata (Argentina) 
for their financial supports. 



[1] K. Mannheim, Astropart. Phys. 3, 295 (1996). 

[2] E. Waxman and J. Bahcall, Phys. Rev. Lett. 78, 2292 (1997). 

[3] IceCube Collaboration Astropart. Phys. 20, 507, (2004). 

[4] P. Jain, S. Kar, D. W. McKay, S. Panda, J. P. Ralston, Phys. ReV. D 66, 065018 (2002). 
[5] L. A. Anchordoqui, C. A. Garcia Canal, H. Goldberg, D. Gomez Dumm, F. Halzen, Phys. 
ReV. D 74, 125021 (2006). 



12 



[6] J. Kwiecinski, A. D. Martin and A. M. Stasto, Phys. ReV. D 59, 093002 (1999). 

[7] N. Arteaga-Romero, C. Carimalo, A. Nicolaidis, O. Panella, G.Tsirigoti, Phys. Letter B 409, 
299 (1997). 

[8] E. Waxman and J. N. Bahcall, Phys. Rev. D 59, 023002 (1999) 

[9] A. Nicolaidis and A. Taramopoulus, Phys. Lett. B 386, 211 (1996). 

[10] Luis A. Anchordoqui, Matthew M. Glenz and Leonard Parker, Phys. Rev.D75, 024011, (2007). 

[11] R. Gandhi, C. Quigg, M. H. Reno, I. Sarcevic, Astropart. Phys. 5 81:110 (1996). 

[12] A. M. Dziewonski and D. L. Anderson, Phys. Earth Planet Inter. 25, 297 (1981) 

[13] E.J.Eichten, K.D.Lane, M.E.Peskin, Phys. Lett. Lett. 50, 811 (1983). 

[14] PDesiati, [astro-ph/0611603| 

[15] M.Ackermann et al, Astroparticle Physics 22, 339-353 (2005). 



13 



